Optimal Topological Test for Degeneracies of Real Hamiltonians 
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We consider adiabatic transport of eigenstates of real Hamiltonians around loops in parameter 
space. It is demonstrated that loops that map to nontrivial loops in the space of eigenbases must 
encircle degeneracies. Examples from Jahn- Teller theory are presented to illustrate the test. We 
show furthermore that the proposed test is optimal. 

PACS numbers: 03.65. Vf, 31.30.Gs, 41.20.Jb, 41.20.Cv 



o 
o 

(N 

Os 

(N 
> 
oo 

o 

O 

m 
o 

a,: 
^— > 

CT 1 - 



X 



Sign reversal of real electronic eigenfunctions when 
continuously transported around a degeneracy was dis- 
covered by Herzberg and Longuet-Higgins hi and was 
later used to construct a topological test for conical in- 
tersections ■ Such intersections are abundant in molec- 
ular systems and are important because they signal 
a breakdown of the Born-Oppenheimer approximation. 
The topological test in has been used to detect coni- 
cal intersections in LiNaK and ozone Q . The sign change 
of the electronic eigenfunctions gives rise to the molec- 
ular Aharonov-Bohm effect 0|, which has recently been 
experimentally tested |(| and theoretically investigated 
0- The structure of adiabatic wave functions and the 
sign change pattern in the vicinity of degeneracies have 
also been analyzed in the context of quantum billiards. 
The behavior of the real wave functions for such systems, 
studied by analog experiments on microwave resonators 
8|, has been interpreted in terms of both the standard 
!j and the off-diagonal E3 geometric phases. The mi- 
crowave resonator experiments have motivated general 
theoretical treatments concerning both the concomitant 
geometric phases and structure of the wave functions 

EES EH. 

In this Letter we study the eigenvectors of a real matrix 
Hamiltonian on a loop in parameter space. We show that 
the behavior of the eigenvectors may imply the presence 
of a degeneracy even if none of the eigenvectors changes 
sign around the loop. This result is a generalization of 
Longuet-Higgins' topological test for intersections 0. It 
is also proven that the generalized test exhausts all topo- 
logical information associated with the behavior of the 
eigenvectors, concerning the presence of degeneracies. 

We begin by the following topological fact. Let X and 
Y be topological spaces. If L is a trivial loop in X, and 
if F : X — > Y is continuous, then F(T) is a trivial loop 
in Y. 

To prove this, note that if G is a homotopy between L 
and a point xq € X, then F o G is a homotopy between 
F(T) and F(x ) G Y. 

Now let H(Q) be annxn matrix Hamiltonian, written 



in the basis {\i}} of the n dimensional Hilbert space H. 
We suppose that H(Q) is real, symmetric, and contin- 
uous for each Q = (Qi, . . . , Qd) in parameter space Q, 
which we assume to be a simply connected subset of R d . 
The eigenvectors of H{Q) can always be chosen real. We 
call the space of real vectors Ureal, which through the 
expansion coefficients in the basis {\i)} can be identified 
with R". The set N rea i of normalized real vectors is the 
sphere S^ 1 in R". 

Consider a simply connected surface S in Q, bounded 
by the loop L. Longuet-Higgins' theorem asserts that 
if a certain eigenvector \^%{Q)) of H(Q) changes sign 
when continuously transported around T, then there is 
a point on S where \ipi(Q)} becomes degenerate with an- 
other state. We note that the theorem implies that if 
H(Q) is nondegenerate on S, then ±\tpi(Q)} represent 
two continuous functions from S to S*™ -1 . 

As a first step towards a generalization of Longuet- 
Higgins' theorem we consider two-level systems. The set 
of normalized vectors N rea i is then the circle S 1 . We 
prove that if an eigenvector of a real two-level system rep- 
resents a nontrivial loop in S* 1 when continuously trans- 
ported along a loop L, then L must encircle a degener- 
acy. The proof is by reductio ad absurdum. Suppose that 
H(Q) is nondegenerate on the surface 5 bounded by L. 
A real eigenvector | ± (Q)) then represents a continuous 
function F from S to N rea i = S . Since L is trivial {S 
being simply connected), so is the loop traced out by 

l±(Q)>- 

As an illustration, let us consider the coupling matrix 
Hamiltonian of the E <£> e Jahn- Teller system 12] 



H(fi,0) 

kp cos 9 
kpsmO 



| gp 2 cos 29 kp sin 9 — | gp 2 sin 28 
lgp 2 sm29 — kpcos9 — \gp 2 cos2# 



(1) 
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where p and 8 are polar coordinates of parameter space 
Q = R 2 , and k and g are the linear and quadratic cou- 
pling strengths, respectively. In this system there are 
four degeneracies: one at the origin p = and three 
at p = 2k/g and 9 = 7t/3, 7t, and 5tt/3. Continu- 
ous transport of an eigenvector around any single de- 
generacy produces a sign change |l2j. Longuet-Higgins' 
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test applied to such a loop would thus imply a de- 
generacy. Consider instead a circular loop T given by 
p ^> 2k jg and 9 e [0, 2n], encircling all four degenera- 
cies. The eigenvectors | — (0)) ss sin0|l) + cos0|2) and 
| + (0)) w — cos0|l) + sin0|2) do not change sign as the 
loop is traversed. However, viewed as elements of the 
circle S 1 , | ± (9)} both make a complete clock- wise turn. 
This means that each eigenvector traces out a nontrivial 
loop in Af re ai — S 1 . In this way the presence of degen- 
eracies encircled by T can be detected. 

For the two-level case we have thus obtained a general- 
ization of Longuet-Higgins' test: also loops along which 
the eigenvectors trace out nontrivial loops in S 1 encircle 
degeneracies on every surface bounded by them. How- 
ever, for n > 3, the space of real normalized vectors 
Nreai = S 71 ^ 1 is simply connected, i.e., it contains only 
trivial loops. Thus, following a single eigenvector along 
the loop is insufficient. To generalize the test in this 
case, we consider instead a complete set of 
normalized eigenvectors of H(Q). 

If H(Q) is nondegenerate on S, then ±\tpi(Q)) are, 
for each i, two continuous and globally defined functions 
from S to Without loss of generality we may as- 

sume that {\ipi(Q))}i=i represents a positively oriented 
orthonormal basis of Tl rea i = R". Every such basis can 
be thought of as an element of the n dimensional rotation 
group SO(n). We may thus define a continuous function 
F : S -> SO(rc) as 

/ (i\MQ)) ■■■ (MMQ)) \ 

F(Q)=\ : •, : . (2) 

\(n\MQ)) ■■■ (n\MQ)) J 

Thus, if H(Q) is nondegenerate on S, then its eigenvec- 
tors represent a continuous function from S to SO(n). 
Analogous to the two-level case it follows that the loop 
in SO(n) traced out by F(Q) as Q varies along T, is a 
trivial loop. By reductio ad absurdum we consequently 
arrive at the following result. If the n eigenvectors of 
H(Q) represent a nontrivial loop in SO(n) when taken 
continuously around T, then there must be at least one 
degeneracy of H(Q) on every simply connected surface S 
bounded by T. 

This result makes it possible to detect the presence of 
a degeneracy by considering eigenvectors on a loop in Q 
even if they do not change sign around the loop. It con- 
stitutes the promised generalization of Longuet-Higgins' 
test. 

We note that SO(2) is homeomorphic to S 1 and is thus 
infinitely connected: it contains one class of nontrivial 
loops for each nonzero integer. The fundamental group 
of SO(2) is the additive group of integers Z. Determin- 
ing whether a loop is trivial amounts to counting the 
number of times the eigenvector encircles the origin. For 
each n > 3, however, SO(n) contains only one class of 
nontrivial loops, that all become trivial when traversed 



twice. The fundamental group of SO(n > 3) is the two- 
element group 7i2 ■ In order to apply the test in this latter 
situation, we need to know how to determine whether a 
loop in SO(n) with n > 3 is trivial or not. In principle, 
this can be done by lifting the loop to the universal cov- 
ering space Spin(n) 0] of SO(n). In practice, though, it 
seems difficult to find a method, that works for all n and 
is easily implemented. Below we show explicitly how to 
use the test in the cases n = 3 and n = 4. We use exam- 
ples from Jahn- Teller theory as physical illustrations. 

First, let us consider the n = 3 case. SO(3) is home- 
omorphic to the closed ball of radius 7r with antipodal 
points on its surface identified 0|. A vector 0v in the 
ball represents a rotation around the unit vector v by the 
angle € [0, 7r]. Antipodal points on the surface must 
be identified since a rotation by the angle 7r is the same 
transformation regardless of whether the rotation axis is 
plus or minus v. A loop in SO (3) can thus be viewed as a 
curve that may exit the closed ball at the boundary, and 
enter again at the antipodal point. The loop is trivial 
if and only if the number of piercings of the boundary 
is divisible by two. Such a piercing is characterized by 
cj> = 7T and that v abruptly changes sign |l5j . 

We apply the method to the linear T ® t? Jahn- Teller 
system, described by the coupling matrix 



H(R) = 




(3) 



where R = (X,Y, Z) parametrizes Q = R 3 . The Hamil- 
tonian is doubly degenerate on eight rays in R 3 that go 
out from the origin and into the middle of each octant. 
At the origin there is a three-fold degeneracy. We con- 
sider a loop r lying completely in the X-Y plane. In 
planar polar coordinates p and 9, the eigenvectors ar- 
ranged according to increasing energy read \ipi(9)) = 
-± ( sin 9\1) + cos 9\2) + 13)) , |<0 2 (0)) = cos 0|1) - sin 0|2), 
and |^3(0)) = -^=(sin0|l)+cos0|2)-|3)). Note that the 
eigenvectors are independent of p and that none of them 
changes sign around any loop that does not pass through 
the origin. The function F along the loop is 



F(9) = 



7T 



sm( 



V 



1 

V2 



cos 9 -A= sin t 
— sin0 -^cos( 

-71 



(4) 



We aim to determine which loops in the X-Y plane that 
map to nontrivial loops in SO (3) under F. The rotation 
angle <fi and rotation vector v of F(9) are given by 



I 

1+ 2 



I 



v(0) = — (-cos 0, sin 0-1,(1 -V^) cos ( 



N 



(5) 



respectively, where A is a normalization constant. Eq. 
(J5J shows that (f)(9) = it only for 9 = n/2. It is also 



FIG. 1: Rotation angle cf> in radians and components of the 
rotation vector v, as functions of 8 also in radians. (f> and v 
are computed from the matrix F(6) of Eq. (JIJ representing 
eigenvectors in the T ® T2 system. The first, second, and 
third component of v are represented by the solid, dotted, 
and dashed curves, respectively. Note that cj> equals n at 
9 — tt/2, and that v changes sign at that point. 



clear that v is continuous everywhere except on the ray 
8 = tt/2, where it abruptly changes sign. A plot of <f> and 
v as functions of 8 appears in Fig. ^ Thus, exactly the 
loops in the X-Y plane that cross the ray 8 — tt/2 an 
odd number of times are mapped to nontrivial loops in 
SO (3). This is a concrete manifestation of the validity 
of the test, since any loop that passes 8 = tt/2 an odd 
number of times must encircle the degenerate subset of 
Q. 

Consider now the case n = 4. Let F(t) = (Fij(t)) 
be a loop in SO(4) parametrized by t S R. The first 
column f(t) = (F 11 (t),F 21 (t),F 31 (t),F il (t)) of F(t) then 
represents an element of S 3 tracing out a loop. Define 
the matrix 



T(f(t)) 

( F n (t) 
-F 21 (t) 
-F 31 (t) 



F 2 i(t) 
Fu(t) 
F 41 {t) 
-F 31 (t) 



F 3 i(t) 
-F 41 (t) 
Fn(t) 
F 2 i(t) 



Fu(t) 
F 3 i(t) 
-F 21 (t) 
Fn(t) 



(6) 



being orthogonal and continuous as a function of / G S 3 . 
It follows that T(t) ee T(f(t)) <G SO(4) since it can 
be continuously connected to the identity matrix, cor- 
responding to / = (1,0,0,0) = §i e S 3 . Furthermore, 
since f(t) is a loop in S* 3 , T(t) is a trivial loop. We con- 
sider the loop in SO(4) given by T(t)F(t). This loop is 
nontrivial exactly when F(t) is, since T(t) is trivial. How- 
ever, the loop T(t)F(t) is simpler than F(t). Carrying 



FIG. 2: The same plot as in Fig. Q but (f> and the vector 
components are calculated from the matrix A(6) of the G®g 
system along the loop V. Since (j> equals tt at one point, V 
must encircle a degeneracy. 



out the matrix multiplication, yields 



T(t)F(t) 



/ 1 \ 

A n (t) A 12 (t) A 13 (t) 

A 21 (t) A 22 (t) A 23 {t) 

\0 ,431 (t) A 32 {t) A 33 {t) J 



(7) 



Thus, the loop T(t)F(t) is trivial if and only if the loop 
A(t) = (Aijit)) in SO (3) is trivial. 

We have applied the above procedure numerically to 
the G £g) g Jahn- Teller system. The coupling matrix for 
this system reads 



H(g) = -qkfV2 

( 93 54 

34 -53 

9i ~ 93 -92 + 94 

\ 92 + 94 gi + 53 



5i - 53 52 + 54 \ 

-52 + 54 5i + 53 
-5i 52 
52 5i / 



(8) 



where g — (qgi, ■ ■ ■ , 954) are normal modes and kg is the 
linear coupling constant. We consider a loop T described 
by 951 = 954 = cos 8 and qg 2 = qg 3 = sin 8, 8 £ [0, 2n]. 
This loop encircles a four-fold degeneracy at the origin 
of parameter space and it can be checked that none of 
the eigenvectors changes sign along V. Fig. [21 shows the 
angle and vector of rotation of the matrix A(8) along 
r. <j) equals tt once, meaning that the loop in SO (4) 
traced out by the eigenvectors is nontrivial. The four- 
fold degeneracy at the origin is thus detected. 

Let us now consider whether there can be a better 
topological test than the present one. It is easy to find 
loops encircling degeneracies that map to trivial loops in 
SO(n). This means that the test does not find every de- 
generacy, so there may be room for improvement. The 
rest of this Letter we devote to an argument showing that 
if we are only allowed to consider eigenvectors on a loop in 
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parameter space Q, and if the only information we have 
about the Hamiltonian is that it is continuous, then there 
is no test that can do better in implying degeneracies. To 
prove this, assume that the eigenvectors {\ipi(Q(t)))}2 =1 
and corresponding eigenvalues {Xi(Q(t))}f =1 are known 
along the loop T parametrized by t e R. Suppose 
that the eigenvectors are nondegenerate along T so that 
Ai(Q(t)) < A 2 (Q(t)) < • • • < X n (Q(t)) holds for all t, and 
that the corresponding loop F(Q(tj) in SO(n) is trivial. 
Under these conditions we show that there always exists a 
Hamiltonian H(Q), continuous and nondegenerate on all 
Q, having exactly the eigenvectors {|V'j( < 9(*)))}?=i an d 
eigenvalues {A;(Q(i))}? =1 on T. Thus, if F{Q{t)) is triv- 
ial, the Hamiltonian can be nondegenerate everywhere, 
meaning that no test can imply a degeneracy. 

For the proof we need to assume that the loop T is 
homeomorphic to S 1 and that there is a homeomorphism 
D from R d onto itself, such that D(T) is the unit circle 
in the X1-X2 plane. (Xi, . . . , Xj) = X = D(Q) denotes 
the coordinates of the image of a point Q € R d under D. 
Such a D exists for all physically interesting Y. 

We consider first the unit disc in the X1-X2 plane, 
parametrized by the usual polar coordinates p G [0, 1] 
and 9 <E [0, 2tt]. By the homeomorphism D we 
may view F(Q(t)) and {Xi(Q(t))}2 =1 as functions of 9, 
i.e., F{6) = F(D- 1 (cose,sm6,0,...,0)) and X,(9) = 
Xi (D _1 (cos 9, sin 6,0, ... ,0)). The eigenvalues can be 
continuously extended to the unit disc preserving the 
nondegeneracy by defining Aj(p, 9) = pXi{9) + (1 — p)Aj, 
where ^4^ are any real numbers satisfying A\ < A 2 < 
...<A n . Let A(p,9) = dmg[Xi(p,9), . . . ,X n (p,9)}. Also 
the eigenvectors can be continuously extended to the disc. 
This is exactly because they represent a trivial loop in 
SO(n). Let the function G : [0, 1] x [0, 1] -» SO(n) be a 
homotopy between the loop F(9) and a constant element 
i?o of SO(n) such that G(0, s) = i? for all s e [0, 1] and 
G(1,9/2tt) = F(9). A continuous extension of F to the 
whole disc is (9) = G(p, 9/2n). 

We are now in a position to construct a suitable Hamil- 
tonian as a function of the coordinate X. Let p take any 
value in [0,+oo). Define 

f F(p,9)A(p,9)F(p,9) T , if p< 1 

tf(X)= „ . (9) 

I F(l,0)A(l,0)F(l,0) r , if p > 1, 

where T denotes matrix transposition. Note that H(X) 
defined in this way is independent of the coordinates 
X 3 , . . . , Xd, that it is continuous and nondegenerate for 
all X, and that it has {\ipi{9))}f =1 as eigenvectors and 
{Xi(9)}2 =1 as eigenvalues on the unit circle in the X1-X2 
plane. To obtain the desired Hamiltonian as a function 
of Q, simply define H{Q) = H(D(Q)). 

In conclusion, we have described a test for degeneracies 



of real Hamiltonians based on the behavior of their eigen- 
vectors on a loop in parameter space. If one considers a 
complete set of eigenstates, one may detect a degener- 
acy even if none of the corresponding vectors changes 
sign. The test works for all real quantum systems with 
finite dimensional Hilbert spaces, including quantum bil- 
liards and Jahn- Teller systems, and could find explicit 
use in quantum chemistry applied to computed eigenvec- 
tors. We show also that no other topological test can do 
better in detecting degeneracies. 

We wish to thank Tobias Ekholm for discussions. The 
work by E.S. was financed by the Swedish Research 
Council. 
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